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Postural body sway displacements for quiet standing sub- 
jects (measured with a new ultrasonic device) are reported. 
Two of the well known strategies for balancing, namely an- 
kle and hip movements were probed. The data is modeled 
using a Fokker-Plank-Langevin stochastic theory. Both ana- 
lytic and computer simulation techniques are employed. The 
Kramers transition rate for a fall is expressed as a function 
of experimental parameters. The root mean square velocity 
is especially important in determining the fall probability. 
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I. INTRODUCTION 

A problem of considerable interest in human body mo- 
tion dynamics is the prediction of the probability to fall 
in a given environment Essential information is 

contained in the small amplitude displacement noise ex- 
hibited by the body during quiet standing. Our purpose 
is to explore how the quiet standing noise data may be 
used to compute the probability (more precisely the tran- 
sition rate per unit time) of falling in various external 
environments. 

In the previous work of other groups, the data on body 
displacement noise were taken using a multicomponent 
force plate upon which the subject quietly stands 
In the present work, we employ a sound wave assessment 
(SWA) device which measures body displacements using 
sound wave echoes; e.g., by measuring the time taken 
for a sound wave pulse to travel from the source to the 
standing subject and/or the time taken for a sound wave 
pulse to travel from the standing subject to a receiver. 
The experimental details concerning the SWA device will 
be discussed in Sec. II. 

Three fundamental movement strategies to maintain 
a standing balance are well known [Q: (i) body move- 
ment from the ankles, (ii) body movement from the hips, 
and (iii) stepping motions. Stepping motions will not 
be considered any further in this work which concerns 
quiet standing. As described in Sec. Ill, we find very dif- 
ferent frequency scales for hip and ankle motions. The 
ankle motions proceed more slowly than the hip motions 
by about a factor of ten. Another feature of the ankle 
motions, also discussed in Sec. Ill, concerns the stability 
angle with respect to the vertical direction beyond which 
a subject, (without bending or stepping) will fall ||. The 



falling instability will be modeled by an appropriate po- 
tential in Sec. IV. Further, the noise itself will be model 
using the Fokker-Planck stochastic equation Q with a 
noise temperature intimately related to the mean square 
velocity fluctuations of the standing subject. 

Apart from the analytic work on stochastic equations, 
we have performed computer simulations of quiet stand- 
ing subjects. In Sec.V we compare the computer simula- 
tions with data taken on subjects using the SWA device. 
The agreement between theory and experiment appears 
satisfactory. 

In Sec. VI we make use of the Kramers potential well 
"escape" formula Q . This determines the transition rate 
per unit time for a subject to fall. The Kramers result 
for the transition rate for a fall is given by 
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where ojq is the frequency of the ankle oscillation mode, 
Q is the quality factor for this mode, a is the critical 
displacement amplitude beyond which the ankle mode is 
unstable, and V is the root mean square velocity which 
is determined by the environmental noise temperature. 

In the discussion Sec. VII, we note that the transition 
rate T k for a fall has an essential singularity Tk — > as 
the root mean square velocity V — ► 0. This implies that 
V is the most sensitive parameter which may be used 
to determine the probability of falling. The implications 
will be discussed in the concluding Sec. VII of this work. 



II. THE SOUND WAVE ASSESSMENT DEVICE 

The SWA device employs two small ultrasonic trans- 
ducers each of which can send or receive ultrasonic sig- 
nals. The first is positioned on a stable laboratory stand. 
The second is attached with belt to the lower back of 
quietly standing subject. Each transducer sends sixteen 
ultrasonic pulses every 3.3 x 10~ 2 sec, which are later de- 
tected by the opposite transducer. The distance between 
the two transducers can then be obtained from the time 
required for the pulses to travel from the sender to the 
receiver. 

A single transducer system has been used in indus- 
try, e.g. by Polaroid for photography applications, in 
order to determine the distance to an object. The sin- 
gle transducer is used to both emit and detect its own 
echo. Our use of two transducers improves the accuracy 
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by having less air absorption as well as less spatial disper- 
sion (fanning out) of the pulses. Our position measure- 
ments are quiet precise; i.e. displacements are measured 
to within an instrumental accuracy of 0.02 cm in a band- 
width of 12 kHz. This yields a displacement noise error 
of Sx « 2 (fim/yHz), Thus, we can record the fine 
movements of a subject (both towards and away from a 
static transducer) as a function of time. 



III. FREQUENCY SCALES FOR HIP AND 
ANKLE STRATEGIES 

Two fundamental body movement strategies Q (to 
maintain a standing balance) are apparent in the SWA 
measurements discussed below. These are the quiet 
standing ankle movements and hip movements, which oc- 
cur in widely different frequency bands. The hip motions 
take place on a time scale of Thip ~ 2 sec, while the ankle 
motions take place on a time scale of T an kie ~ 30 sec. 
In order to measure ankle motions, the durations of our 
measurements were rather long, i.e. 3 min, while the 
shorter runs of 1 min duration were sufficient to study 
hip motions. The short times scales, of less than 20 sec 
were sufficient to verify the fractal diffusion estimates of 
hip motions which have appeared in the pioneering stud- 
ies P,§-0 using the force plate technology. 



forth characteristic of ankle motions along with the more 
rapid (and in this case) much smaller hip motions. 

Were the quiet standing subject to move without hip 
motions, i.e. move rigidly, then the stiff body would be 
held steady, if and only if the angle d formed with the 
vertical were less then some critical angle d c . The geo- 
metrical cone d < d c is called the sway envelop or the 
cone of stability ||. Typical values for the critical (for- 
ward) angle is d c « 8° |ll]]. For angles d > $ c , the 
ankles can no longer hold the subject upright. In the ab- 
sence of a step or other supports (such as arm motions) 
a fall will take place. If one considers the displacement 
of the part of the body at a height h above ground, then 
a = htawdc represents the critical displacement ampli- 
tude beyond which the ankle mode is unstable. Such crit- 
ical displacements are normally measured using standard 
tests of the subjects forward or backward reach Jl2|,[l3| . 



IV. FOKKER-PLANCK-LANGEVIN THEORY 

The Fokker-Planck-Langevin stochastic approach || to 
quiet standing swa y p rocesses has been shown to be ex- 
tremely useful JlJ,]l^|. In what follows we apply this 
method to the ankle movements. As discussed in Sec. Ill, 
the ankle movements can be described by a slow oscilla- 
tory motion modulated by faster random noise. 

To mathematically describe the cone of stability, we 
model the metastable potential by the potential 
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FIG. 1. SWA measurement of x(t) for a quiet standing sub- 
ject. The large broad peaks and valleys are from ankle mo- 
tions, and the smaller more rapid oscillations are from hip 
motions. 



To illustrate the decomposition of movements into an- 
kle motions and hip motions we plot, in FIG.l, the dis- 
placement as a function of time for a typical quiet stand- 
ing subject. One notes the slowly varying drifts back and 
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where ujo is the frequency of the ankle oscillation mode, 
m is the mass of the mode described by the displacement 
x, and a is the critical displacement. 

The ankle mode, as shown in the experimental FIG.l, 
is far from a perfect oscillation. There exists a random 
force f{t) and a finite quality factor Q for the mode 
which, in the Langevin theory, are intimately related. 
The equation of motion for the mode may be written as 



W \ . 2 
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+ U'(x) = f(t)+F(t), (3) 



where F(t) is an applied force, and f(t) is a "white noise" 
random force obeying 

</(*)/(0>^=(^^)*(*-t'), (4) 

and where T n is noise temperature. The noise tempera- 
ture is by no means equal to the ambient temperature of 
the room in which the subject stands. Rather the noise 
temperature describes all of those fluctuations which cou- 
ple into the coordinate x. In particular, all of the envi- 
ronmental and internal fluctuations which couple into x 
contribute to the root mean square velocity V 
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which can be obtained from experimental data such as 
that pictured in FIG.l. The noise temperature T n is here 
defined in terms of the root mean squared velocity V via 
the equipartition theorem fllq] 
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Thus, for the case in which the applied force is written 
F(t) = m(3(t), the random force is written f(t) = ma(t) 
and the potential U(x) — m4>(x), the Langevin equation 
reads 

x + (j^j x + wlx + tfix) =a(t) +/3(t), (7) 
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Eqs.(7), (8) and (9), with a hip modulation given by 

P(t) = PmaxCOs(ui hlp t) (10) 

can be used to perform a computer simulation of the 
ankle mode of motion. 



For the healthy quiet standing subject, shown in 
FIG.l, an output of the random motion simulation is 
shown in FIG. 2. The qualitative similarity between a 
number of quiet standing subjects and simulations allow 
us to conclude that the Langevin theory of postural sway 
is reasonable. 

The Langevin theory can also be analytically expressed 
as an equation for the probability P(x,v, t)dxdv for the 
subject to have a velocity in the interval dv and a dis- 
placement in the interval dx. For example, without the 
hip modulation, the Fokker-Planck equation reads 
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The Fokker-Planck and the Langevin formulations of the 
problem are equivalent. The former is useful for ana- 
lytical calculations while the later is useful for computer 
simulations. 



V. POSTURAL SWAY DATA 

Not all of the quiet standing subjects measured ex- 
hibited large undamped ankle mode oscillations. For 
some quiet standing subjects the ankle movements were 
strongly damped. Shown in FIG. 3 is an example of a 
subject with suppressed ankle movements. The hip os- 
cillations are still clearly visible. 



3 
2 



o 

x, 
-1 



■ h \ 








\ A 




- t\! 








v \ 
\ ) 




/ 


\ f \ 


/ \ 


V W 




V 1 






) f 











50 100 150 

(t/sec) 



147 




50 100 150 

(t/sec) 



FIG. 2. Simulated quiet standing posture sway with a fre- 
quency luq = 0.125 sec -1 , a critical reach of a = 12 cm, a noise 
temperature corresponding to V = 0.22 cm/sec, a quality fac- 
tor of Q = 20, a hip modulation frequency of uihi P = L5 sec -1 , 
and a hip modulation amplitude of f3 max = 0.37 cm/sec 2 . 



FIG. 3. Shown is a quiet standing subject with a strongly 
damped ankle mode. The hip oscillations about the diffusive 
ankle motions are clearly visible. 
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The computer simulations for the overdamped case 
proceed exactly as previously described. However the 
quality factor of the mode is taken to be much smaller 
than for the underdamped case. For the quiet standing 
subject, shown in FIG. 3, an output of the random motion 
simulation is shown in FIG. 4. We again conclude, for the 
overdamped case, that the Langevin theory of postural 
sway is reasonable. 
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FIG. 4. Simulated overdamped quiet standing posture 
sway with a frequency luq = 0.125 sec -1 , a critical reach 
of a = 12 cm, a noise temperature corresponding to 
V = 0.31 cm/sec, a quality factor of Q = 0.1, a hip mod- 
ulation frequency of uJhi V = 1.5 sec -1 , and a hip modulation 
amplitude of fimax = 0.56 cm/ sec 2 . 
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FIG. 5. SWA measurement of ^Aa;(r) 2 ^ for a quiet stand- 
ing subject taken from the raw data in FlG.l. The undamped 
oscillation of the ankle mode is clearly visible. The hip mo- 
tions are smoothed out in the mean square displacement fluc- 
tuation. 

In FIG. 6 we have plotted (Aa;(r) 2 ) for the simulation 
in FIG. 2. We note that the oscillations for underdamped 
ankle movements are clearly present although the ampli- 
tude of the oscillation arc high when compared with the 
experimental FIG. 5. However, we still conclude that the 
Langevin model is qualitatively reasonable. 



The ankle movement is more clearly seen in the mean 
square displacement fluctuation defined in mathematical 
Brownian motion theory as 

<Ax(r) 2 ) = 



lim 
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\x(t + r) -x(t)\ dt. 
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In physical systems the time T of a run is finite. In our 
case T = 180 sec. We can then plot (Ax(t) 2 ) on the 
interval < t < 140 sec. 

For the underdamped experimental data of FIG.l, we 
show in FIG. 5 a plot of (Aa;(r) 2 ). The ankle movements 
show an oscillation in a very clear form when the averag- 
ing procedure of Eq.(12) is performed. The hip move- 
ment modulations are largely smoothed away by this 
same averaging procedure. 

The data produced by the computer simulations can 
also be averaged according to Eq.(12) and then compared 
with mean square fluctuations taken from the experimen- 
tal data. The hip movement modulations are barely no- 
ticeable in (As(t) 2 }. 
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FIG. 6. Simulated mean square fluctuation (Ax(t) 2 } with 
parameters the same as in FIG. 2. 



For the overdamped case, the data in FIG. 3 gives rise 
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to an experimental (As(r) 2 ) shown in FIG. 7. Notice for 
the overdamped case the absence of oscillations in the an- 
kle motion. The behavior of (Ax(t) 2 ) on the time scale 
shown is qualitatively not very far from simple Brownian 
diffusion. 



4 



3 - 

CXI 

E 
u 

A 

CM 2 - 




1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — ' 
50 100 

(x/sec) 



FIG. 7. SWA measurement of (Ai(t) 2 ^ for a quiet stand- 
ing subject taken from the raw data in FIG. 3. The oscillation 
of the ankle mode overdamped. 

This diffusion-like behavior is also present in the com- 
puter simulation shown in FIG. 8. The simulated ampli- 
tudes are again large compared with the experimental 
amplitudes. 
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FIG. 8. Simulated mean square fluctuation (Ax(t) 2 } with 
parameters the same as in FIG. 4. 



The Fokker-Plank-Langevin theory is mainly of quali- 
tative significance when comparing simulations to exper- 
imental data. 

VI. KRAMERS "ESCAPE RATE" FOR A FALL 

Although the form of the metastablc potential in 
Eq.(2) is commonly used, the potential has not yet been 
shown to be unique. The qualitative features required 
are a minimum near x = and a barrier to falling near 
x = a. Eq.(2) is plotted in FIG. 9. 
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FIG. 9. The model potential U(x) for describing a fall. 

A quiet standing subject exhibits oscillations about the 
potential minimum. The random force at a noise tem- 
perature T n will at some time push the displacement to 
values | a; | > a over a potential maxima at which time 
there is a fall. The potential barrier protecting the sub- 
ject from a fall is given by 

AU = U(\a\)-U(0) = (^f^j. (13) 

Neglecting hip motion modulations, the Kramers equa- 
tion (jj for the escape rate contains the Boltzman factor 
for overcoming a barrier; It is 

r»-(^) .-/«■-. (14, 

In terms of experimental parameters 

We used the parameters in Eq.(15) to describe twenty 
healthy subjects (between sixteen and forty five years of 
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age) in street clothing with eyes wide open. The reach 
test length scale a was calculated from the critical angle 
i? c given in the literature . We find the Kramers pre- 
dicted "times for a fall" to be anywhere from one week to 
three years. The time shortens considerably as the mean 
square velocity increases. Of course, this prediction by 
no means implies that these subjects will all experience a 
fall within the next three years! The fall transition rate 
of our model will be further discussed in the next and 
concluding section. 

VII. DISCUSSION 

We have exhibited measurements of displacement noise 
with quiet standing subjects along with computer simu- 
lations based on a Fokker-Plank-Langevin model. For 
all ranges of the quality factor Q, from underdamped 
Q » 1 to overdamped Q << 1, the Langevin equa- 
tion computer simulations were in qualitative agreement 
with the data. In this regard, one should note that the 
same subject when measured on different days can ex- 
hibit both overdamped and/or underdamped behavior, 
as well a somewhat different effective noise temperatures 
T„, i.e. a somewhat different root mean square velocity 
V. 

Since the Kramers transition rate for a fall depends 
sensitively on the root mean square velocity V, it follows 
that a persons transition rate for a fall can vary from 
day to day depending on the "noise temperature" . For 
example, on a day when a person is tired we conjecture 
that a fall is more likely than on a day when the person 
is wide awake. However, we presently have no data with 
which to prove this conjecture. 

While the above Fokker-Planck-Langevin theory seems 
to be qualitatively correct, the theory overestimates the 
displacement noise amplitudes and thereby also overesti- 
mates the transition rate for a fall. Apart from leaving 
out the step strategy for avoiding a fall, we believe that 
we have perhaps also underestimated the subtlety of the 
hip strategy. Thus far, the hip motions have been de- 
scribed by an added modulation force in our computer 
simulations, and this inclusion gives at least qualitative 
agreement with experimental data. Nevertheless, the hip 
modulations might be more strongly correlated with the 
center of mass body coordinate |lj],[l8| than is presently 
being included. Further improvements on the present 
model are presently being pursued on both the experi- 
mental and theoretical level. 
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